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The Pierre Auger Observatory detects ultra-high energy cosmic rays by implementing two complementary air- 
shower techniques. The combination of a large ground array and fluorescence detectors, known as the hybnd 
concept, means that a rich variety of measurements can be made on a single shower, providing much improved 
information over what is possible with either detector alone. In this paper I describe the hybrid reconstruction 
approach and the latest hybrid measurements. 



1. Introduction 

The energy spectrum of primary cosmic rays 
extends over more than 12 decades in energy. It 
follows, over a large range, a simple power law 
indicating its non-thermal nature. The highest 
energy cosmic rays (which means particles with 
energies above 10 EeV = 10^^ eV) must be pro- 
duced by the most extreme non-thermal process 
in the Universe. Their energy spectrum, arrival 
directions, and composition can be inferred from 
air shower observations. However, the origin of 
these particles is still unknown and there is not 
a consensus whether there is a steepening in the 
energy spectrum at around 60 EeV. Determining 
the shape of the spectrum in this energy range 
in sufficient detail requires an enormous detec- 
tor exposure. The Pierre Auger Observatory was 
designed to study these ultra-high energy cosmic 
rays. Although under construction, the Auger 
Southern site is the largest cosmic ray detector in 
the world. 

A fundamental characteristics of the Pierre 
Auger Observatory is its capability of hybrid re- 
construction of cosmic ray showers Qj. Two in- 
dependent detectors, the Surface Detector (SD), 
which samples the shower particles at ground, 
and the Fluorescence Detector (FD), which col- 
lects the fluorescence light emitted by the shower 
particles along their path in the atmosphere, are 
able to measure the energy and direction of the 
same cosmic ray shower. The enhanced capabili- 



ties of the Auger hybrid detector are examined in 
this paper. 

The Pierre Auger Observatory was designed to 
observe, in coincidence, the shower particles at 
ground and the associated fluorescence light gen- 
erated in the atmosphere. This is achieved with 
a large array of water Cherenkov detectors cou- 
pled with air-fluorescence detectors that overlook 
the surface array. It is not simply a dual exper- 
iment. Apart from important cross-checks and 
measurement redundancy, the two techniques see 
air showers in complementary ways. A single air 
shower is detected 3-dimensionally. The ground 
array measures the 2-dimensional lateral struc- 
ture of the shower at ground level, with some 
ability to separate the electromagnetic and muon 
components. The fluorescence detector records 
the longitudinal profile of the shower during its 
development through the atmosphere. 

A hybrid event is an air shower that is simul- 
taneously detected by the fluorescence detector 
and the ground array. The Observatory was orig- 
inally designed and is currently being built with 
a cross-triggering capability. Data are recovered 
from both detectors whenever either system is 
triggered. If an air shower independently triggers 
both detectors the event is tagged accordingly. 
There are also cases where the fluorescence detec- 
tor, having a lower energy threshold, promotes a 
sub-threshold array trigger. Surface stations are 
then matched by timing and location. This is an 
important capability because these sub-threshold 
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hybrid events would not have triggered the array 
otherwise. The geometrical reconstruction of the 
air shower's axis is accomplished by minimizing a 
^ function involving data from all triggered ele- 
ments in the eye and on the ground. The recon- 
struction accuracy is far better than the ground 
array counters and the single eye could achieve 
independently [2]. 

The combination of the air fluorescence mea- 
surement and the particle detection on the ground 
provides an absolute energy calibration. The FD 
measurements determine the longitudinal devel- 
opment profile, whose integral is proportional to 
the total energy of the electromagnetic particle 
cascade. The SD independently estimates the 
shower energy by evaluating the particle density 
at 1000 m from the core. This ground parame- 
ter, called S'(IOOO), is estimated from the surface 
stations by fitting the measured densities to a lat- 
eral distribution function 3 . A conversion factor 
that relates 5(1000) to the shower energy (based 
on the FD information) is extracted from hybrid 
events. This reduces significantly the dependence 
on air shower models and on assumptions of the 
primary composition. 

The Pierre Auger Collaboration started the 
construction of the Southern site in 2002 in the 
city of Malargiie, which is located at an eleva- 
tion of 1400 m a.s.l. in the province of Mendoza, 
Argentina. After a successful operation of a pro- 
totype experiment 0], the Southern Observatory 
started operation in hybrid production mode in 
January, 2004. Surface stations have a 100% duty 
cycle, while fluorescence eyes can only operate on 
clear moonless nights. Both surface and fluores- 
cence detectors have been running simultaneously 
14% of the time. The number of hybrid events 
represents 10% the statistics of the surface ar- 
ray data. The Southern site will be completed in 
2007. Another important objective is to obtain a 
uniform exposure over the full sky. This will be 
achieved by constructing a second instrument in 
the Northern hemisphere. Each site is conceived 
to cover an area of 3000 km^ in order to collect 
about 1 event per week and site above 60 EeV. 
The Northern Observatory is planned to be sited 
in the U.S. 



2. Hybrid Geometrical Reconstruction 

A hybrid detector achieves the best geometri- 
cal accuracy by using timing information from all 
the detector components, both FD pixels and SD 
counters. Each element records a pulse of light 
from which one can determine the central time of 
the pulse and its uncertainty. Each trial geome- 
try for the shower axis yields a prediction for the 
times at each detector component. Differences 
between actual and predicted times are weighted 
using their corresponding uncertainties, squared, 
and summed to construct a value. The hy- 
pothesis with the minimum value of is the re- 
constructed shower axis. 
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Figure 1 . Light track of a hybrid event as seen by 
two adjacent fluorescence telescopes. The full line 
is the fitted shower-detector plane. The squares 
represent the surface detectors that also triggered 
in this event. 



In the FD, cosmic ray showers are detected as 
a sequence of triggered pixels in the camera. An 
example of an event propagating through two ad- 
jacent FD telescopes is presented in Fig. The 
first step in the analysis is the determination of 
the shower-detector plane (SDP). The SDP is the 
plane that includes the location of the eye and the 
line of the shower axis. (See the sketch in Fig. [21) 
Experimentally, it is the plane through the eye 
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which most nearly contains the pointing direc- 
tions of the FD pixels centered on the shower axis. 
(See fitted line in Fig.^) Using a known axis pro- 
vided from the Central Laser Facility (CLF), de- 
scribed in Ref. 0, the SDP reconstruction error 
can be evaluated by comparing the space angle 
between the normal vector to the experimentally 
determined SDP and the known true normal vec- 
tor. This uncertainty in the SDP is not greater 
than about 0.1°. 




Figure 2. Illustration of the geometrical shower 
reconstruction from the observables of the fluo- 
rescence detector. 



Next, the timing information of the pixels is 
used for reconstructing the shower axis within the 
SDP. As illustrated in Fig. [3 the shower axis can 
be characterized by two parameters: the perpen- 
dicular distance Rp from the eye to the track, and 
the angle ■0 that the track makes with the hori- 
zontal line in the SDP. Each pixel which observes 
the track has a pointing direction which makes an 



angle Xi with the horizontal line. Let to be the 
time when the shower front on the axis passes the 
point of closest approach Rp to the eye. The light 
arrives at the ith pixel at the time 

R 

U=to + ^cot[{i; + X^)/2]. (1) 
c 

The shower parameters are then determined by 
fitting the data points to this functional form. 
Using the FADC electronics, such a monocular 
reconstruction may achieve excellent accuracy. 
However, the accuracy of the monocular recon- 
struction is limited when the measured angular 
speed dx/dt does not change much over the ob- 
served track length. An example is shown in 
Fig. 121 For these events (usually short tracks) 
there is a small curvature in the functional form 
of Eq. ^ such that there is a family of possible 
(Rpjip) axis solutions. This is also illustrated in 
Fig. 21 Not only the uncertainties in Rp and ip 
are large, but also they are strongly correlated. 
This translates directly into an uncertainty in 
the other shower parameters, especially in the 
reconstructed shower energy. This asymmetric 
uncertainty in the energy and angular resolution 
are important drawbacks of the monocular recon- 
struction. 

The aforementioned degeneracy can be broken 
by combining the timing information from the SD 
stations with that of the FD telescopes. This is 
called the hybrid reconstruction. The hybrid so- 
lution for the example shown in Fig. |21 is shown 
in Fig. ^ as a white star and its uncertainty as 
the smaller (full) ellipse. 

Since the SD operates with a 100% duty cy- 
cle, most of the events observed by the FD are 
in fact hybrid events. There are also cases where 
the fluorescence detector, having a lower energy 
threshold, promotes a sub-threshold array trig- 
ger. Surface stations are then matched by tim- 
ing and location. This is an important capability 
because these sub-threshold hybrid events would 
not have triggered the array otherwise. In fact, 
the time of arrival at a single counter at ground 
can suffice for the hybrid reconstruction. 

The reconstruction uncertainties are evaluated 
using events with known geometries, i.e. laser 
beams. Since the location of the CLF (approx- 
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Figure 3. Functional form that correlates the 
time of arrival of the light at each pixel with the 
angle between the pointing direction of that par- 
ticular pixel and the horizontal line within the 
shower-detector plane. 



imately equidistant from the first three fluores- 
cence sites) and the direction of the laser beam 
are known to an accuracy better than the ex- 
pected angular resolution of the fluorescence de- 
tector, laser shots from the CLF can be used to 
measure the accuracy of the geometrical recon- 
struction. Furthermore, the laser beam is split 
and part of the laser light is sent through an opti- 
cal fiber to a nearby ground array station. Thus, 
the axis of the laser light can be reconstructed 
both in monocular mode and in the single-tank 
hybrid mode. The resolution of the monocular 
and hybrid reconstructions are compared in Fig.[3 
for the distance between the eye and the CLF, 
and in Fig. for the angle of the axis. The re- 
sults are very encouraging. With the monocular 
reconstruction, the location of the CLF can be 
determined with a resolution of ^ 500 m. After 
including the timing information of the single wa- 
ter tank, the resolution improves by one order of 



Figure 4. Solutions (and la regions) for the axis 
fit shown in Fig. |3| The large uncertainty (and 
strong correlation) of the monocular reconstruc- 
tion is broken using the timing information from 
the surface detectors (Hybrid). The stars mark 
the solutions that minimize the for the axis 
reconstruction. (Following the notation of this 
paper, xo = 180° - -0-) 



magnitude with no systematic shift. 

As mentioned before, the laser light from the 
CLF produces simultaneous triggers in both the 
surface and (three) fluorescence detectors. The 
recorded event times are used to measure and 
monitor the relative timing between the two de- 
tectors. The time offset between the first fluo- 
rescence eye and the surface detector is shown in 
Fig. This time offset is crucial for the accuracy 
of the hybrid reconstruction, and it has been mea- 
sured to better than 50 ns [HI. The contribution 
to the systematic uncertainty in the core location 
due to the uncertainty in the time synchroniza- 
tion is 20 m. 

Using the timing information from the eye's 
pixels together with the surface stations to re- 
construct real air showers, a core location reso- 
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Figure 5. Difference between the reconstructed 
and true distance from the eye to the vertical laser 
beam using the monocular and hybrid techniques. 
The location of the laser is known to 5 m. 



Figure 6. Angular difference between recon- 
structed and true direction of the laser beam us- 
ing the monocular and hybrid techniques. The 
laser beam is vertical within 0.01°. 



lution of 50 m is achieved. The resolution for 
the arrival direction of cosmic rays is 0.6° 
These results for the hybrid accuracy are in good 
agreement with estimations using analytic argu- 
ments [7|, measurements on real data using a 
bootstrap method and previous simulation 
studies in]. 

3. Energy Determination 

After successful reconstruction of the event ge- 
ometry, the FADC traces of the FD pixels are 
analyzed in order to obtain the light emitted 
along the shower axis. An atmospheric scatter- 
ing model is used in this step to transform the 
light received at the FD back to the light emitted 
from the shower axis. The geometrical height, as 
observed by the telescopes, is converted to gram- 
mage of atmosphere. The amount of fluorescence 
light emitted from a volume of air is proportional 
to the energy dissipated by the shower particles in 
that volume. The observed longitudinal light pro- 
file represents the energy loss in the atmosphere, 
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Figure 7. Time offset between the surface detec- 
tor and one of the fluorescence detectors. The 
variation in time is within the uncertainty of 
50 ns. 



which in turn is proportional to the number of 
charged particles in a given volume. The result 
for a hybrid shower of 6* ~ 60° is shown in Fig. |S1 
(bottom). The line represents the best fit to a 
Gaisser-Hillas function yielding a primary 

energy oi E = (23 ± 6) EeV, in good agreement 
with the S'(IOOO) determination. 
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Figure 8. An air-shower event detected both in 
water tanks (top: array view), and in a fluores- 
cence detector (bottom: longitudinal profile). 



4. Observatory Status 

The status of the Southern Observatory on 
May 10, 2006 is summarized in Fig.|5| There were 
953 fully operational tanks at that time. The first 
two fluorescence sites (Los Leones and Coihueco) 
were fully operational, i.e. running six telescopes 
each, in June, 2004. The third site (Los Morados) 
started operation on March 18, 2005. The fourth 
and last site (Loma Amarilla) is currently under 
construction and is scheduled to start operation 
by the end of 2006. The present average rate is 
50 hybrid events per night per eye, for a total of 
~ 40000 events up to June, 2006. At this rate, 
4000 hybrid events per month are expected when 
the Southern site is completed. 

5. Hybrid Measurements 

In the Pierre Auger Observatory two com- 
plementary experimental approaches are com- 
bined on a shower-by-shower basis within one sin- 
gle experiment. Such redundancy allows cross- 
correlations between experimental techniques. 



Figure 9. Status of the Pierre Auger Southern 
Observatory on May 10, 2006. There were 1138 
deployed tanks and 953 of them were fully oper- 
ational. Also shown are the locations of the flu- 
orescence detectors with their fields of view (for 
the three fully operational detectors) . The fourth 
fluorescence site (Loma Amarilla) is under con- 
struction. 



thereby controlling the systematic uncertainties. 
Furthermore, there is an improvement in the res- 
olution of the energy, mass, and arrival direction 
of the reconstructed primary particles. Data are 
being used to verify both the performance of the 
individual components |llll2j . as well as to pro- 
duce the highest quality reconstructed air-shower 
events Results are very promising |[14j and 

underline the advantages of the hybrid approach. 

5.1. Energy Spectrum 

Hybrid air shower measurements are utilized 
in the spectrum analysis to avoid dependence on 
specific numerical simulations of air showers and 
detector responses to them. The analysis is also 
free of assumptions about the primary nuclear 
masses. The FD provides a nearly calorimetric, 
model-independent energy measurement: fluores- 
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cence light is produced in proportion to energy 
dissipation by a shower in the atmosphere. Hy- 
brid data establish the relation of shower energy 
to a ground parameter. Moreover, hybrid data 
determine the trigger probability for individual 
tanks as a function of core distance and energy, 
from which it is found that the SD event trig- 
ger is fully efficient above 3 EeV for zenith angles 
less than 60°. The SD exposure is then calculated 
simply by integrating the geometric aperture over 
time. It is the continuously operating surface ar- 
ray which provides the high statistics with unam- 
biguous exposure. Assigning energies to the SD 
event set is a two-step process. The first step is to 
assign an energy parameter, 5*38, to each event. 
Then the hybrid events are used to establish the 
rule for converting 5*38 to energy. 

The energy parameter 5*38 for each shower 
comes from its experimentally measured 5'(1000). 
It may be regarded as the 5(1000) measurement 
the shower would have produced if it had arrived 
38° from the zenith. As it can be seen in Fig. 1101 
S'as is well correlated with the FD energy mea- 
surements of high quality hybrid events. The fit- 
ted line gives an empirical rule for assigning en- 
ergies (in EeV) based on S^s (in VEM): 

E ^O.Wx Slii°^. (2) 

The SD acceptance is not saturated below 3 EeV. 
But the hybrid events used in Fig.^| which start 
at ~ 1 EeV, are those with core locations and ar- 
rival directions such that they have a probability 
greater than 0.9 for satisfying the SD trigger and 
quality conditions. The distribution over energy 
produced by this two-step procedure becomes the 
cosmic ray spectrum presented in Ref. |15| . 

5.2. Anisotropy 

Due to the much improved angular accuracy, 
the hybrid data sample is ideal for anisotropy 
studies and, in particular, for point source 
searches. Results on a search for a point -like 
source in the direction of the Galactic Center 
(GC) using these hybrid events were presented in 
Ref. ■ These events have a better angular reso- 
lution (0.7° at 68% c.l. in the EeV energy range). 
Considering the hybrid events with an energy be- 
tween 10^^'^ eV and 10^^'^ eV, no significant ex- 




Log(S38/VEM) 



Figure 10. Shower energy (measured from 
the longitudinal profile) as a function of the 
ground parameter S^s- Each point is a hybrid 
event recorded when there were contemporaneous 
aerosol measurements, whose longitudinal profile 
included shower maximum in a measured range of 
at least 350 g cm~^, and in which there was less 
than 10% Cherenkov contamination. The fitted 
line is log(£;) = -0.79 + 1.06 log(5'38). 



cess is seen in the GC direction. For instance, in 
an optimal top-hat window of 1.59cr ~ 0.75° ra- 
dius, 0.3 events are expected while no single event 
direction falls within that circle. This leads to a 
source fiux upper-bound at 95% c.l. of 

$f = ^ 0.15 km-2 yr-^ (3) 

(The factor ^ is close to unity and parametrises 
the uncertainties in the flux normalization.) 

Note that the energy assignments of the FD 
apply regardless of the assumed primary compo- 
sition (except for a small correction to account for 
the missing energy) , be they protons or heavy nu- 
clei. However, the acceptance has a dependence 
on composition because different primaries de- 
velop at different depths in the atmosphere. Since 
a quality requirement for hybrid events is to have 
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the maximum of the shower development inside 
the field of view of the telescopes, this affects the 
sensitivity to different primaries. The bound ob- 
tained is indeed conservative even if the bulk of 
the cosmic rays are heavy nuclei. 

5.3. Composition 

An upper limit for the photon fraction in cos- 
mic rays with energies above 10^^ eV was derived 
using the hybrid data set ^7]. The importance 
of the photon fraction lies in the fact that top- 
down models predict a considerable proportion 
of photons among the generated particles. The 
measured photon flux is thus a valuable indicator 
of these non-accelerator models. 

The method used in Ref. (I7j to distinguish be- 
tween hadrons and photons in the hybrid data 
exploits the information on the longitudinal pro- 
file of the air shower. In fact, this is the first 
such limit on photons obtained by observing the 
fluorescence light profile of air showers. The at- 
mospheric depth at the shower maximum, Xmax, 
is commonly used as a discriminant observable for 
the cosmic ray composition because lighter nuclei 
penetrate, on average, more deeply into the at- 
mosphere. Above 10^^ eV, showers initiated by 
photons develop significantly deeper in the atmo- 
sphere than hadronic showers. 

No candidate for a primary photon was found 
in the hybrid data taken between Jan, 2004 and 
Feb, 2006. By comparing the observed X^ax of 
each hybrid event to predictions from hadronic 
simulations, an upper hmit of 16% (at 95% c.l.) 
was derived. This results confirms and improves 
the existing limits above 10 EeV. 

This analysis is currently limited mainly by the 
the small number of events. The number of hy- 
brid events will considerably increase in the next 
few years, and much lower primary photon frac- 
tions can be tested. An upper limit of ~ 5% could 
be achieved, for example, with two more years 
of data taking. Moreover, the larger statistics 
will allow to increase the threshold energy above 
10 EeV where even larger photon fractions are 
predicted by the models. A similar limit 15%) 
but at higher energy (e.g. above 40 EeV), would 
be well below existing limits and severely con- 
strain non-acceleration models. 



6. Conclusions 

The Pierre Auger Observatory is a hybrid de- 
tector with excellent capabilities for studying the 
highest energy cosmic rays. Much of its capabil- 
ity stems from the accurate geometric reconstruc- 
tion it achieves. The shower geometry is recon- 
structed combining information from the eyes and 
the ground detectors. Arrival directions are de- 
termined to a small fraction of a degree and the 
shower core is located to an accuracy of about 
60 m [nE|. 

The construction of the Southern Observatory 
is well under way. Eighteen FD telescopes and 
more than 60% of the surface array are in oper- 
ation taking data routinely. At the present rate 
of deployment, construction will be finish in mid 
2007. Detectors are performing very well and the 
first results are very encouraging. 

Emphasis is placed on hybrid analysis that pro- 
vide unprecedented quality in geometry, energy, 
and mass reconstruction. Of utmost importance 
for the near future will be the determination of 
the energy spectrum to study the GZK feature, 
and the search for anisotropics in arrival direc- 
tions. It is important to note that both SD and 
FD techniques have different systematics, and re- 
sults are preliminary at this stage while the Ob- 
servatory is under construction. The possibility of 
studying the same set of air showers with two in- 
dependent methods is valuable in understanding 
the strengths and limitations of each technique. 
The hybrid analysis benefits from the calorimetry 
of the fluorescence technique and the uniformity 
of the surface detector aperture. 

In parallel to the completion of the Southern 
Observatory and to the analysis of data towards 
the first scientific results, R&D has started for the 
development of the Northern site in the US. 
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